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Abstract 

We address the issue of edge detection in Synthetic Aperture Radar imagery. In 
particular, we propose nonparametric methods for edge detection, and numeri- 
cally compare them to an alternative method that has been recently proposed 
in the literature. Our results show that some of the proposed methods display 
superior results and are computationally simpler than the existing method. An 
application to real (not simulated) data is presented and discussed. 

Keywords: Edge detection, hypothesis testing, image analysis, multiplicative 
noise, small samples, speckle, SAR 



1. Introduction 

Synthetic Aperture Radar (SAR) images are an important source informa- 
tion in many applications, such as urban planning, environmental monitoring, 
crop management, oil prospection, mining exploration, wind detection, animal 
life detection, among others. A SAR is a coherent radar of high resolution that 
works on-board using a synthetic antenna of a movable platform, like an airplane 
or a satellite, covering extended surfaces and producing images. SAR systems 
employ the Doppler effect and processes the signal obtaining high spatial reso- 
lution in the direction of the platform motion [24] . 

During the data collection, the target remains illuminated under the antenna 
beam for a few moments and is observed by the radar from positions induced 
by its movement throughout the platform trajectory. The radar illuminates the 
target with a succession of pulses of a given frequency. The energy is propagated 
in all directions, and part of it returns to the antenna (this return is called 
'echo'). The sensor measures both the intensity and the delay between the 
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signals sent from and received back by the antenna. The image is then formed 
based on the energy returned by each point on the surface. 

Some of the most important features of the SAR sensor for remote sensing 

are: 

• The active nature of the instrument makes it independent of other iUumi- 
nation sources, being able to gather images at any time. 

• Microwaves intensity are not significantly affected by the presence of clouds, 
so image acquisition is possible in most metheorological conditions and in 
regions with permanent cloud coverage. 

• SAR images can have high spatial resolution, e.g. of less than one meter, 
thus making the study of small scale phenomena possible. 

• SAR images contain complementary information to that provided by op- 
tical images. The selection of frequency band, polarization and angle of 
incidence in SAR imagery allows the discrimination of different surface 
properties. 

A SAR sensor emits and receives electromagnetic waves of complex nature 
and, therefore, the received signal can be stored in different formats: complex, 
intensity, amplitude and phase [22]. 

In a SAR image, it is possible to distinguish several types of roughness or 
texture, according to which one can classify the different types of covers: 

• Homogenous areas: Surfaces of very little texture; for example, crops, 
deforestation, and, under some conditions, snow, water or ice. 

• Heterogeneous areas: Surfaces that display some texture; for example, 
forests on not very pronounced reliefs, among others. 

• Extremely heterogeneous areas: Surfaces with intense texture; for in- 
stance, urban areas, among others. 

"Texture" , in the context of SAR imagery, should be understood as a mea- 
sure of the number of objects in a cell of the size the of the wavelength employed 
by the sensor. A fine texture corresponds to a large number of objects per cell, 
while coarse or extremely heterogeneous textures are those for which only a few 
objects are counted per cell. The Japanese Earth Resources Satellite JERS-1, 
for instance, operates on L-Band (1.3 GHz, 23.5 cm wavelength) and the Euro- 
pean Remote Sensing Satellites ERS-1 and ERS-2 use C-Band (5.3 GHz, 5.6 cm 
wavelength) . 

A problem of paramount importance in the analysis of images is segmenta- 
tion: the process that divides an image in its constituent parts or objects. Its 
main goal is to group image areas that have similar characteristics. One of the 
basic principles in the segmentation process is the detection of discontinuities. 
Edges are the borders of the objects and are therefore quite useful for their 
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segmentation, registration and identification. Edges can be thought as the lo- 
cations where abrupt changes in intensity or in other important characteristic 
occur. 

The quahty of SAR images is degraded by speckle, a degradation which 
follows from the use of coherent illumination, i.e., when the signal phase is em- 
ployed in the image formation. Such degradation is characteristic of technologies 
that employ microwaves, sonar, laser and ultrasound. 

The presence of speckle makes edge detection difficult, since most algorithms 
identify regions using local characteristics. Though speckle should not be re- 
garded as noise, since it has a deterministic nature and is reproducible, from 
the image practitioner viewpoint it can be considered a random effect and can 
be conveniently described by stochastic laws; c.f. Richards (2J, Sec. 4.3.1]. It is 
not convenient to only use pointwise information when detecting edges under 
speckle; it is necessary to analyze the image using sets of pixels that provide 



local information [15 1. 



Different approaches can be used to locate the edges between regions in a 
SAR image. A particularly attractive and well performing statistical method 
was proposed by Gambini et al. It is based on the family of Q distribu- 

tions, which can be successfuly used to describe areas with different degrees of 
homogeneity [HI, [H, [H, [ll] . 

In this work we consider intensity imagery, described by the Gj law. This 
distribution is indexed by the number of looks L > 1, the scale parameter 
7 > and the roughness parameter a < 0. The former can be controlled when 
generating the image or in postprocessing stages, and is a measure of the signal- 
to-noise ratio [2J, Sec. 4.3.1]. The value of the roughness parameter is of interest 



in many applications, since it can be used as an indicator of land type. The 
scale parameter relates to the relative power between the reflected and incident 
signals fll|. 

Figure [1] presents three different targets and the corresponding values (or 
range of values) of the roughness parameter (a). Small values of a (e.g., 
a < — 10) are associated with homogeneous areas, such as pastures. Values of 
a e [—10, —4] are characteristic of heterogeneous regions, for example forests. 
Finally, larger values of a (say, — 4 < a < 0) are observed in extremely hetero- 
geneous areas, such as urban regions (gI. [l9|. 
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Our chielf goal is to develop and assess new methods for edge detection 
in SAR images using computationally efficient nonparametric statistical infer- 
ence. Gambini et al. [3| showed that a method based on maximum likelihood, 
which is presented in Section [3.11 is more precise than four commonly used tech- 
niques: two based on raw data (maximum discontinuity and fractal dimension) 
and two based on estimates (maximum discontinuity and anisotropic smoothed 
roughness). Our numerical evaluation reveals that some of the edge detection 
strategies we propose in this paper outperform the best method available to 
date, i.e., maximum likelihood at considerably lower computational cost. 

The paper unfolds as follows. Section [5] presents the model. Section [3] 
discusses the use of nonparametric statistics for SAR image detection, while 
Section |4] presents the results. An application that uses real (not simulated) 
data is presented and discussed in Section [5] Conclusions and directions for 
future work are outlined in Section [51 



2. The multiplicative model 

The physics of SAR image formation leads to the multiplicative model: the 
random field Z describing the observations can be viewed as the product of two 
independent random fields, which are not observed directly, namely: X and Y . 
The former {X) models the properties of the imaged area (backscatter), whereas 
the latter {Y) models the speckle noise due to the use of coherent illumination. 

Speckle noise in intensity L-looks format follows a Gamma distribution, de- 
noted by F ~ r(L, L), whose density is given by 

Jviy) = Y^D^"^'^ expi-Ly}, L > 1, y > 0. 

In SAR images modeling, the smallest value of L is 1, which corresponds 
data with the highest spatial resolution, at the expense of lower signal-to-noise 
ratio. The number L can be assumed to be known or estimated beforehand 
from the entire image. We shall assume that the number of looks is known. 

The backscatter exhibits different degrees of homogeneity, and different mod- 
els have been used to encompass this characteristic. The reciprocal of Gamma 
distribution is a useful and tractable model P, 0, 11 El, [H, [3, El, H [H, [IS 



2l|,|26j. The random variable X follows this law, denoted AT ~ F ^(a,7), if its 



density function is 

fx[x)^—^ ra;~"~^exp|--|, -a,7,a;>0. 

7"! \—a) L X ) 

Let X ~ r~^(a,7) and Y ~ T{L,L) be two independent random variables. 
It can be shown that the random variable Z — XY follows a distribution, 
denoted Z ~ ^j(a,7,i), whose density is 



7"r(L)r(-a) i-/ + Lzy 
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where — a > is the roughness parameter, 7 > is the scale parameter and 
L > 1 is the number of looks • 

In this work we consider independent samples. Correlated fields are treated 
by Bustos et al. @. 

The r-th noncentral moment of Z is given by 



E[Z''] 



T)Y[L- 



r(-a)r(L) 



(2) 



if —a > r, and 00 otherwise. 

The distribution is very attractive for modeling data with speckle noise, 
due to its mathematical tractability and because it is able to describe informa- 
tion from most types of areas. 

Figure m shows 5j(a,7a,3,3) densities, where ja,L is the value of 7 that 
delivers unit expected value for given a < — 1 and L. The densities are presented 
in semilogarithmic scale, showing that they have heavy (linear) tails with respect 
to the Gaussian distribution which displays quadratic behavior. It is noticeable 
that larger values of a lead to larger variances; in fact, the variance is not finite 
when a > — 1. 




Figure 2: G^{q, 7^,3, 3) densities for a = —3, —6, —12 (dots, dashes, solid) 

Figure 13] shows nine patches of speckled data for one, three and eight looks. 
Each patch consists of 3 x 3 images with varying roughness and contrast: each 
column shows images with same mean (1, 5 and 10, left to right) while rows 
from bottom to top show images with a G {—2, —5, —20}. It can be readily seen 
that finding edges in speckled data can be a hard task due to the existing local 
variation. 
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(a) Single look 



(b) Three looks 



(c) Eight looks 



Figure 3: Patches of speckled data 



Gambini et al. [14 , 15 1 used an analogy estimator based on moments of order 
1/2 and 1. Since we wish to extend their work, we shah use the same approach. 

Let (Zi, . . . , Zn) be a vector of independent identically distributed random 
variables, with common distribution ^x(a, 7, i), a < —1/2, 7 > and L known. 
Define the r-th sample moment as 



1 " 
V = - V 



n 

i=l 



with r = 1/2 and 1 (that is the reason why we assume a < — 1 rather than 
a < 0). 

From equation ([2]) it is immediate that 

fi\ r(-a-i)r(L + i) 



and that 



Replacing the population moments by their sample counterparts, and the 
parameters by the corresponding estimators, we arrive at the following system 
of two equations: 

/7\ r(-S-i)r(L + i) ^ ^ 
= r(-a)r(L) ' 

and 

-7\^ r(-5-i)r(£ + i) 
^^J r(-S)r(L) ' 

which leads to the following equation that can be solved numerically in order 
to obtain an estimator for a: 

miT{-a)T{L)L mlV^ {-a)T'^ {L) L 



r(-a-i)r(L + i) r2(-s- i)r2(L + 1)' 
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By plugging the value of a into equation ^ we obtain 7. 



Frery et al. 1^ showed that computing ML estimators for the family is 



prone to severe numerical instabilities, and they proposed an iterative algorithm 
that alleviates this problem. Recently Pianto and Cribari-Neto [2^ analyzed 
this issue and found that it is related to a flattening of the likelihood function, 
and they proposed a correction based on resampling. Cribari-Neto et al. Q 
showed that ML estimators for that distribution can be quite biased, and they 
evaluated the effectiveness of improving them by several resampling techniques. 



Vasconcellos et al. [26| proposed an analytical bias correction for ML estima- 



tors; they showed that there is a wide range of practical situations for which the 
corrected estimator effectively reduces both bias and mean square error of the 
original ML estimator. Bustos et al. @ derived M-estimators and showed that 
such robust estimators are superior to the classical estimators in the presence 
of corner reflectors, a common source of contamination in SAR images. Al- 
lende et al. [1| derived AM-estimators (M estimators with asymmetric influence 
functions), motivated by the shape of the Ga density. Overall, their estimators 
outperform both ML and M-estimators. 

In the next section we shall summarize the main techniques that are available 
for edge detection in SAR imagery, with special emphasis on those that explicitly 
employ statistical models and techniques. 



3. Edge detection in SAR imagery 

Statistical edge detection is described in Bovik et al. These authors 
introduced nonparametric statistics for edge detection under Gaussian additive 
noise. They showed the usefulness of the median and the Wilcoxon-Mann- 
Whitney tests for edge detection with the help of a sample image interpreted 
visually. 

Fesharaki and Hellestrand proposed an algorithm for edge detection by 
using a t-test, while Beauchemin et al. Q used a nonparametric alternative 
based on the Wilcoxon-Mann- Whitney statistics for detecting changes between 
adjacent pixel neighborhoods. Although these tests may be appropriate for 
specific types of images, they may not detect changes in local grey level values 
in images with low signal-to^oise ratios, as is the case of SAR imagery. 

Hoon Lim and Ju Jan g Il6l| compared two-sample tests for edge detection in 
noisy images. Later Lim [l7| described a new edge detector based on the robust 
rank-order test, an alternative to the Wilcoxon test, using £ x £ windows for 
detecting all possible edges in noisy images. This method is based on testing 
whether an ^ x ^ window is partitioned into two sub-regions. 

The detection of edges in images with speckle noise has been studied by 
many authors, and a variety of techniques are presented and compared by 
Gambini et al. 15 1, whose proposed approach outperforms all competing 
procedures with an acceptable computational cost. Their method is based on 
fitting contours of objects and regions using B-splines; such curves depend on a 
few parameters and can be easily computed from control points and smoothness 
conditions 
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Gambini et al. 14, 1^ used amplitude data for finding edges between dif- 
ferent regions. Such data, if squared, follow the 5j distribution. Under this 
model, the regions of the image with different degrees of homogeneity are char- 
acterized by the parameters of the distribution. If a point belongs to the edge 
of the object, then a sudden change in the parameter values is expected in its 
neighborhood. 

In what follows we present five techniques for detec ting edges in SAR images, 



namely the method proposed by Gambini et al. IJ, |15| , and four alternative 
techniques based on ranks. The two main advantages of nonpar ametric tests 
are that 

1. they do not depend upon the data distribution, and 

2. they are robust to extreme observations. 

These features are of particular interest in the case at hand, namely, speckled 
imagery. 



3.1. The Gambini algorithm for edge detection 

This is an iterative procedure that refines an initial region until the final 
result is achieved. It is based on the fact that if a point belongs to the object 
edge, then a sample taken from its neighbourhood should exhibit a change in 
the parameter values and, therefore, should be considered a transition point. 

Consider N image segments, s''^\i £ {1,...,N}, of the form s*^*^ = CPi, 
where C is the centroid of the initial region, the extreme Pi is a point outside 
the region and 9i = Z(s(*\ s^'"*"^)) is the angle between two successive segments, 
as shown in Figure 21 




Figure 4: Radial lines from the centroid C to the exterior of the region, with a separation of 

Consider a strip of pixels around each segment s*^*), as illustrated in Figure[5l 
partitioned into two areas by a candidate edge point Ci (the small red dot). 
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The observations within the strip are taken to come from two different models, 
namely one corresponding to the object C/§(a^, 7^, L) (the dark area in Figure[5]) 
and the other from the background Q^(ar,^r, L). 




Figure 5: Segment, strip and candidate edge point (small red dot) 

The parameters {ae,ji) and (a^iTr) index the region and its background, 
respectively, and their estimation is described in Section [21 note that the esti- 
mators depend upon the the transition point q. 

In order to find the transition point on each segment s*-*-*, an objective func- 
tion is considered: the sample likelihood, which is given by 

j m 

^(af,7^,ar,7r) = ]^Pr(zj;a£,7£) X J| Pr(2;i; a^, 7r), 

i—l 

where j is the number of observations in the strip around segment s^'^ lying 
between C and q, and m is the number of observations in the strip lying between 
Ci and Pi. 

In order to find the transition point, we maximize the log-likelihood function 

j m 

C = \n{e) = ^\nfgo{zf,ai,-fi)+ ^ In /g;n (z,; a^, 7r) 

i=l i=j + l 
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for every possible value of j. Using equation ([1]), and assuming that L is known, 

7rT(L)r(-S?)(f, + Lz,)« 

7r"'^r(L)r(-s;)(> + iz,)i-^ 

Finally, the estimated transition point on the segment is given by 

J — argmax£. (5) 

i 

Figure IHl shows typical values of the objective function, taken along a straight 
line segment. 



Figure 6: Values of the objective function for a segment of straight 



3.2. Nonparametric edge detection 

In the typical two-sample situation, the experimenter has two samples from 
possibly different populations, and wishes to use a statistical test to determine 
whether the null hypothesis that the two populations are identical should be 
rejected. That is, the experimenter wishes to identify differences between the 
two populations on the basis of their random samples. 

An intuitive approach to the two-sample problem is to combine both samples 
into a single ordered sample, and assign ranks to the observations, regardless 
the source population. A possible test statistic is the sum of the ranks assigned 
to each population. If such a sum is small (large), there is indication that the 
values from that population tend to be smaller (larger) than the values obtained 
from the other population. The null hypothesis should be rejected if the sum of 
the ranks associated with one sample is considerably larger than the other sum. 
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3.2.1. The Mann- Whitney test 

The data consist of samples from two populations. Let Xi, X2, . . . , X„ de- 
note the random variables sampled from population A and Yi, I2, • • • , Ym denote 
the variatcs obtained from population B. Assign the ranks 1, . . . , N = n + m 
to the individual observations in the combined sample. Let R{Xi) and R{Yj) 
denote the ranks assigned to Xi and Yj for all i and j. If two or more sample 
values arc equal (ties), assign to each the average of the ranks that would have 
been assigned had there been no ties. 

The following assumptions are made: 

• Both samples are random samples from the respective populations. 

• In addition to independence within each sample, there is mutual indepen- 
dence between the two samples. 

• The measurement scale is at least ordinal. 

Let F{x) and G{x) be the distribution functions of X and Y, respectively. 
We wish to test 



In many situations, as in edge detection, differences between distributions 
imply that Pt{X < Y) does not equal 1/2. Therefore, we can rewrite the above 
hypotheses as 



The null hypothesis can be tested using the Mann- Whitney test, which is 
unbiased and consistent. The test statistic is computed as follows. When there 
are no or just a few tics, the sum of the ranks assigned to the sample from 
population A can be used as a test statistic: 



When the number of ties is large, one can subtract the (null) mean from T and 
divide the resulting difference by the (null) standard deviation to get 



■Ho : F{x) = G{x) 
■Hi : F{x) G{x). 



■Ho : Pr(X <Y) = - 
Ui : Pr(X <Y)^\. 



n 



T = ^i?(X,). 
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3.2.2. The Kruskal-Wallis test 

The Mann- Whitney test for two independent samples was extended to the 
problem of analyzing k independent samples, fc > 2, by Kruskal and Wallis in 
1952. The experimental situation is that k random samples have been obtained 
from k possibly different populations, and one wishes to test the null hypothesis 
that all populations are identical. The Kruskal-Wallis test statistic is a function 
of the ranks of the observations in the combined sample (like the Mann- Whitney 
test statistic). 

The data consist of k random samples of possibly different sizes. Denote 
the ith random sample of size rij by Xi\,Xi2, . . . , . Then, the data may be 
arranged in columns: 

Sample 1 Sample 2 • • • Sample k 

Xl^2 X2,2 Xk^2 



Xl^m ^2,712 Xk,nk 

Let N denote the total number of observations, i.e., N = X^iLi""-!- Assign 
rank 1 to the smallest of the N observations, rank 2 to the second smallest 
observation, and so on. Let Ri be the sum of the ranks assigned to the ith 
sample: 

Ri = Y,Ri^ij) i = l,2,...,fc. (6) 

Compute Ri for each sample. Under ties, assign the average rank to each tied 
observation, as in the previous test. 
The following assumptions are made: 

• All samples are random samples from their respective populations. 

• In addition to independence within each sample, there is mutual indepen- 
dence amongst samples. 

• The measurement scale is at least ordinal. 

• Either the k population distribution functions are identical or else some of 
the populations tend to yield larger values than the remaining populations. 

The null and alternative hypotheses are 

Hq : All k population distribution functions are identical. 
Hi : At least one population tends to yield observations larger than 
at least one of the remaining populations. 

Since the Kruskal-Wallis test is designed to be sensitive against differences 
among means in the k populations, the alternative hypothesis is sometimes 
stated as 
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7^1 : The k populations do not all have identical means. 
The test statistic T is defined as 

r^^l>:--i^y (T) 




where 




Under no ties, S'^ simplifies to N{N + 1)/12, and the test statistic reduces 

to 

1 9 ^ f>2 

N(N + 1) ^ rii ^ ' ^ ' 

When the number of ties is small or moderate, there is little difference between 
equations ([7]) and ([5]), and equation ([5]) is preferred. 

3.2.3. The squared ranks test for variances 

The squared ranks test can be used to assess equality of variances across two 
or more independent, random samples which have been measured using a scale 
that is at least interval [?']. 

The data consist of the two random samples. Let Xi, . ■ , denote 
a random sample of size n from population A and Yi,l2j • • • ,Ym represent a 
random sample of the size m from population B. Define 

= \Xi - fill, i = 1, . . . ,n, 

and 

Vj = 1^3 -^■2], j = l,...,m, 

where /Ii and /X2 are the two sample means. 

Assign ranks 1 to n + m to the combined sample, as usual. If several values 
of U and/or V are equal (ties), assign to each the average of the ranks that 
would have been assigned to them had there been no ties. 

The required assumptions can be stated as follows: 

• Both samples are random samples from the respective populations. 

• In addition to independence within each samples, there is mutual inde- 
pendence between samples. 

• The measurement scale is at least interval. 
The null and alternative hypotheses are 

Hq : X and yare identically distributed, except for possibly different means, 
■Hi : Var(X) ^ Var(F). 
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If there are no values of U tied with values of V, the sum of the squared 
ranks assigned to population A can be used as a test statistic: 



V[i?((7,)] 



If there are ties, subtract the (null) mean from T and divide the difference by 
the (null) standard deviation to get 

T-nW 



nm v-^iv T-,4 njn f p2\2 

N{N-1) l^i=\ N-iy^ I 

where N ^ n + m, and i?^ is the average of the squared ranks of both samples 
combined: 

-J I n m 



also, 



N n m 

^^[i?(c/0]^ + 5]TO)]4. 

i=l i=l j=l 

3.2.4- The TPE empirical statistic 

The empirical statistic TPE is also based on ranges, and is well suited for 
situations where one wishes to test whether two samples come from the same 
distribution or from distributions with different means and/or variances. The 
data consist of the two random samples. Let Xi,X2, . . . , Xn denote a random 
sample of size n from population A and let Yi , I2 , ■ • ■ , Ym denote a random 
sample of the size m from population B. Assign ranks 1, . . . , N = n + m, and 
compute X = n"^E^=l^W^: Y = m-^T,f=i(Y)„ and De = \X -Y\. The 
mean rank of the combined sample is — {N + l)/2. The empirical statistic 
TPE is £: = \De - 

3.3. Proposal 

Our chief goal is to perform edge detection in SAR images using the afore- 
mentioned nonparametric tests statistics instead of the likelihood function pre- 
sented in equation 

The new noparametric edge estimates j on segment s*-*-* are given by 

• Mann- Whitney estimate: J = arg maxj Ti j- , 

• Kruskal-Wallis estimate: argmaxj T^j-, 

• Squared Ranks estimate: J= argmaxj Tyj, and 

• TPE empirical estimate: '] = argmin^ Ej, 
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where j denotes the edge candidate coordinate. 

Figure [7] shows values of the of the Mann- Whitney, Kruskal-Wahis, Squared 
Ranks and the TPE Empirical test statistics taken along a straight line segment 
of simulated data. The corresponding positions of the maxima (Mann- Whitney, 
Kruskal, Variance) or minimum (TPE) are taken to be the transition point 
between the two regions. Data generation was carried in such a way that the 
true transition point is at j — 50. We note that all four methods successfuly 
identify the edge point. 




20 40 60 80 100 20 40 60 80 100 

j j 




20 40 60 80 20 40 60 80 100 



Figure 7: Test statistic values: Mann- Whitney, Kruskal- Wallis, T. Variances and TPE 



4. Numerical results 

This section presents simulations performed to obtain a quantitative assess- 
ment of the performances of the proposed edge detection methods. All sim- 
ulations were run on personal computers with Intel© Pentium© IV CPUs of 
3.20 GHz running Windows XP operating system. The programming language 
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used was Ox version 4.10 [for details, seeQ. All graphics were produced using 
R version 2.6.1 



We shall evaluate the errors made when estimating the edge point (local 
error) in several parametric situations. For each situation, 1,000 simulated 
rectangular windows of sizes 20 x 100 (20 rows, 100 columns) are filled with 
with samples from Qj distributions. 

Each window is composed of two halves, and we consider all possible com- 
binations of roughness parameters —ae 6 {3,4,6,8,10,12,14,16,18,20} and 
-ar e {2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20}, and number of 
looks L G {1,3,8}. Here, ag (ar) denotes the value of a used for data gen- 
eration to the left (right) of the true edge point. These values span a variety 
of images often encountered in practice, ranging from single-loook to smoothed 
multi-look. The scale parameter 7 equals 

T{-a)T{L)L 
^"'^ " T{-a-l)r{L + l)' 

so that the distribution mean equals one in all situations. 

It is noteworthy that edge detection algorithms are therefore evaluated in a 
very difficult situation, namely, when the areas have the same mean and only 
differ in texture. 

Situations in which ai — ar are not considered, since they entail no edge. 
We thus consider ((10 x 19) — 10) x 3 = 540 {ag, ar and L) cases. For each of 
them, we simulate 1, 000 windows. In each of these 540, 000 windows, the edge 
is detected by the five techniques already described, and the error is defined as 
the absolute difference between the true edge (which is located at position 50) 
and the detected edge. 

The errors are stored in an 1000-dimensional array DmU) defined in the 
following way: 

DMU) = \50-PTij)\, J = 1,2,..., 1000, 

where Px{j), j — 1, . . . , 1000, is the transition point identified by method AI in 
the jth sample. Note that M ranges in the set 

{Gambini, Kruskal, Mann- Whitney, Variance, TPE}. 

Let f{M) denote the percentage of times in which the estimated edge is more 
than 5 pixels apart from the true value, using method Af , i.e., 

^ #{je{l,...,1000}:iJM(j)>5} 
•'^ ' 1000 

These are the error rates we report. Method Mi is considered to be more accu- 
rate than method Mj whenever f{Mi) < f{Mj), i.e., whenever the percentage 
of errors of method Mi is smaller than that of method Mj , except for errors up 
to 5 pixels. 

Figure [8] depicts four of the situations assessed: three where single (L = 1) 
look data are contrasted (Figures 8(a)[|8(b7] and 8(c) ), and one with L ^ 8 looks 
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(Figure [8(d)| ). The single look data show the differences between heterogeneous 
(ag — —8, Figure 8(a) I, and homogeneous {ae = —12, Figure [8(b)] ) ; ai = —18, 



Figure [8(c)] ) areas and strips of varying roughness. The L = 8 data only presents 
the difference between a homogenous area {ae = —18) and several other strips 
(Figure [8(d) I . The areas to the left are formed by independent draws from the 
Qxicte, 7af ,L, L) distribution. The strips to the right are formed by independent 
outcomes of the t/j(ar, 7q^,l, i) law with = {—4, —6, —8, —10, —12, —20}. 
The contrast among regions has been enhanced in order to make visualization 
easier; actual data are harder to differentiate. 




Figure 8: Four situations assessed, varying roughness to the left and strips ar = 
{-4, -6, -8, -10, -12, -20} to the right 



Figure [S] summarizes the main results of our study, regarding error rates. 
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Each column presents the error rates as a function of L, the number of looks, 
provided a roughness for the left half of the strip ae- We notice that the error 
rates are consistent with respect to L, since their behavior does not alter sig- 
nificantly. Different lines, for a given L, show the behavior of the error rates as 
a function of the roughness. Mann-Whitney is the only nonparametric tech- 
nique that exhibits poor performance, but only when ai > ar] all other proce- 
dures are competitive. Kruskal edge detection behaves similarly to Gambini 
procedure, with minor differences that do not consistently benefit either. TPE 
and Variance are slightly worse than Kruskal and Gambini, mainly when 
multilook {L = 3, 8) imagery is used. The error rates are consistently larger 
when the left and right roughnesses are similar; notice that the curves peaks 
shift to the left from top to bottom. Overall, the Kruskal and Gambini error 
rates decrease as the number of looks increases, the same pattern holding for 
the degree of heterogeneity. 

Figure[TU]presents the average execution times for the ai — —3 and L = 1 sit- 



uation, which is representative of all remaining execution times. As Figure 10(a) 



shows, the time spent by Gambini decreases as the value of ar increases, but it 
is consistently larger than the nonparametric techniques detection times. The 
differences in execution times are of three orders of magnitude and occur in all of 
the situations herein assessed. Figure [TO(b)| shows the average execution times 
of the nonparametric methods, which never exceeded one second; Kruskal was 
the fastest followed by Mann- Whitney. 

As noted earlier, the hardest task one can face is the discrimination of regions 
indexed by similar parameter values. Table [T] presents a subset of our numerical 
results and focuses on the most challenging cases. It is clear from the figures in 
this table that Kruskal is consistently faster than Gambini, reaching a scale 
of 10^ in some cases. The error rates of the two methods are similar. 



5. Application to real data 

Figure [Tl] presents a SAR image for which the estimated number of looks 
is 3.2. The window considered has 101 x 181 pixels, and was obtained over 
agricultural fields of Oberpfaffenhofen, Germany, in the L-band, by the ESAR 
sensor jSH]. 

The original data displayed almost no differences, and it was enhanced for 
visual purposes. Two main regions can be observed, namely: the dark (left) 
and light (right) areas. The estimated roughness parameter in both areas are 
approximately equal to 7.5 (i.e., a — 7.5), thus implying that both areas are 
slightly heterogeneous, probably due to relief. 

An edge was detected in each of five non-overlapping strips of 20 x 181 
pixels. The Mann- Whitney, Kruskal- Wallis and Gambini estimates agree and 
are indicated by red dots in the middle of each strip. The estimates, thus, 
provide accurate starting points for any subsequent edge detection algorithm. 
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Table 1: Error rates and execution times of challenging situations 
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6. Conclusions 

Our chief goal was to propose alternative techniques for edge detection in 
speckled i mag ery. The proposed methods were compared to that of Gambini 
et al. [3, Il5l | The techniques here assessed do not try to eliminate existing 
speckle, but to extract information from its statistical properties. The model 
used to describe these data is the distribution, which Mejail et al. 3 3 
show can be used as an universal model. 

The methods under assessment aim at identifying edges between regions 
with different degrees of roughness, which, in turn, is determined by a, the 
roughness parameter of the Gxio^i li distribution. Homogenous, e.g., pastures, 
heterogeneous regions, e.g., forest, and very heterogeneous regions, e.g., urban, 
targets are considered. 

In order to compare the performances of the five methods, a Monte Carlo 
experiment was carried out. Two criteria were employed in the comparison, 
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namely: the error made in detecting an edge and the execution time. The 
thoughest setting was considered, namely, when the areas on both sides of the 
edge have the same mean and number of looks, differing only on the roughness. 

From the experimental results it was observed that the Kruskal detector 
performed slightly better than Gambini. Compared with the TPE, Variance 
and Mann- Whitney detectors, the Kruskal performed much better. 

It is noteworthy that all detection methods perform well when the edge 
separates areas with very different degress of roughness. 

When the generated image represents homogenous zones, the error frequen- 
cies of all methods are relatively high, nevertheless this problem is alleviated by 
increasing the number of looks L, causing the edge detection to become more 
accurate. 

The Kruskal detector displays the best results, both with respect to error 
and execution time. The latter is, in many cases, 1000 times smaller than that 
of Gambini. 

The TPE method for edge detection performs well when ae differs from 
by at least three units. 

Mann- Whitney has good performance (both in terms of error rate and 
speed) when ai < ar- When > ar the method is unable to locate the edge, 
mostly because it only accounts for the ranks of one of the samples. 

Overall, our results suggest that the Gambini edge detection technique can 
be successfully replaced by our Kruskal. By doing so, one achieves similar 
precision at a much lower computational cost (the latter is approximately one 
thousand times faster). 

A promising line of research is the use of stochastic distances, as derived by 
Nascimento et al. |21|. 
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(a) ae = -8, L = 1 



(b) ae = -8, L = 3 



(c) ae = -8, L = 8 




(d) ai = -12, L = 1 



(e) at = -12, L = 3 



(f) at = -12, L : 



(g) 



-18, L = 1 



(h) at : 



-18, L = 3 



(i) at = -18, L : 



Figure 9: Error rates as functions of the number of looks {L £ {1, 3, 8}, left to right columns) 
and of the roughness (ai G {—8, —12, —18}, top to bottom rows) 
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(a) All techniques (b) Nonparametric techniques 

Figure 10: Average execution times (in seconds) for on = —3 with L = 1 



Figure 11: ESAR image and detected edge points 
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